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We discuss a geometrical method to define a preferred reference frame for precessing binary 
systems and the gravitational waves they emit. This minimal-rotation frame is aligned with the 
angular-momentum axis and fixes the rotation about that axis up to a constant angle, resulting in an 
essentially invariant frame. Gravitational waveforms decomposed in this frame are similarly invariant 
under rotations of the inert ial frame and exhibit relatively smoothly varying phase. By contrast, 
earlier prescriptions for radiation-aligned frames induce extraneous features in the gravitational- wave 
phase which depend on the orientation of the inert ial frame, leading to fluctuations in the frequency 
that may compound to many gravitational- wave cycles. We explore a simplified description of post- 
Newtonian approximations for precessing systems using the minimal-rotation frame, and describe 
the construction of analytical/numerical hybrid waveforms for such systems. 



I. INTRODUCTION 

One of the central goals of modern numerical relativ- 
ity is the accurate simulation of compact binary systems, 
in particular the computation of the gravitational wave- 
forms emitted by these systems. These waveforms pro- 
vide crucial input into the construction of accurate tem- 
plate banks necessary for detection and parameter esti- 
mation based on matched filtering [1, 2] in gravitational- 
wave detectors such as LIGO, Virgo, LCGT [3-6], and 
possible space-based detectors such as LISA [7-9]. More 
generally, detailed and accurate knowledge of waveforms 
provides a dictionary to relate measured waveforms to 
the specific details of the astrophysical sources that give 
rise to those waveforms, allowing gravitational- wave ex- 
periments to fulfill their proper role as tools for extremely 
high-precision astrophysics. 

For such a bank of gravitational waveforms to be use- 
ful, however, it must not be restricted to an astrophys- 
ically unrealistic subset of the space of source param- 
eters. Numerical-relativity simulations must eventually 
treat binary systems with a broad range of mass ratios, 
spin magnitudes, and spin orientations. Many of the 
fundamental challenges on the first two points have now 
been overcome.^ However, precessing binaries remain a 
formidable challenge. The parameter space of precessing 
binaries is vastly larger than for non-precessing binaries, 
and its exploration is just getting underway. Further- 
more, the absence of precession allows simplifying as- 
sumptions about the properties of the gravitational wave- 
forms, greatly easing post-Newtonian comparisons and 
simplifying gravitational- wave data- analysis strategies. 

When binary systems of black holes (or neutron stars 
or other compact objects) have spin angular momenta 



See reviews [10, 11], as well as [12-14] for simulations that push 
the limits of large mass ratios and spins. 



that are misaligned (with one another or with the orbital 
angular momentum of the pair), the plane of orbital mo- 
tion inclines and precesses over time. In post-Newtonian 
theory, this phenomenon can be interpreted as arising 
from spin-orbit couplings. In fully nonlinear general rel- 
ativity, these dynamical effects cannot easily be under- 
stood in gauge-unambiguous language. However, the ef- 
fects of this precession can be seen unambiguously in the 
waveform, as modulations that trade energy content be- 
tween the various spherical-harmonic modes [15-17]. 

These modulations present difficulties for cataloging 
the gravitational waveforms. In non-precessing simula- 
tions, the standard practice has been to decompose each 
spherical harmonic component of the waveform into a 
time- varying amplitude and phase. Both of these ele- 
ments, in non-precessing cases, can be accurately approx- 
imated before or after merger as simple polynomials or 
exponentially damped polynomials, greatly simplifying 
their description. In precessing systems, this is no longer 
the case. 

Aside from the complexity of fitting the precessing 
waveforms, there is also the concern that a precessing 
system does not have a preferred frame of inertial coor- 
dinates. An overall rotation of the inertial coordinates 
would transform the various waveform modes, modulat- 
ing them in different ways. This means that comparisons 
between waveforms must account for a rotation between 
the inertial frames in which they are measured. In non- 
precessing cases, this difficulty is avoided by the existence 
of a preferred, fixed axis (the axis along which the radi- 
ation is preferentially beamed). This axis is intuitively 
associated with the normal to the orbital plane. 

Schmidt et al. [18] pointed out that this preferred axis, 
while no longer fixed in precessing cases, nonetheless 
still exists and can be used to rotate the inertial spa- 
tial coordinates over time, demodulating the waveform. 
They suggested a method for finding that axis, as did 
O'Shaughnessy et al. [19] more recently. Essentially, the 
idea is to find a rotation operator at each instant in time 
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to maximize the z component of the angular momentum 
in the radiation.^ The z axis of the original frame, ro- 
tated by this operator, is taken to be the preferred axis. 
We refer to this axis as the radiation axis. 

We expect this insight to be important for understand- 
ing and cataloging generic gravitational waveforms. Defi- 
nition of the radiation axis is the first step toward frames 
adapted to precessing binaries, in which the gravitational 
waveforms have simple structure. However, fixing the 
axis does not fix the frame, because of the ambiguity 
in rotations about the axis. From a more formal per- 
spective, the (directed) axis being tracked lives in a two- 
dimensional space: the space of unit vectors in R^, which 
is topologically the two sphere S'^. However the space of 
available rotations, the group manifold of 50(3), is topo- 
logically RP^, which is three dimensional. In general, no 
mathematically preferred method exists to infer a unique 
path in RP^ from a path in S'^; additional conditions 
must be imposed. 

Therefore, the second step toward adapted frames is to 
fix this rotation about the radiation axis. References [18] 
and [19] address the first step by providing suitable def- 
initions of the radiation axis, but deal with the second 
step only implicitly through the choice of parameteri- 
zation of the rotation matrices. A rotation around the 
radiation axis changes the phase of each gravitational- 
radiation mode by an integer multiple of the rotation 
angle. An unsuitable choice of this angle will induce un- 
physical variations in the gravitational-wave phases in 
the adapted frame, even for vanishingly small precession, 
as demonstrated in Sec. IV. 

The present paper addresses the question of how to fix 
the rotation about the radiation axis. Our construction 
is geometric, and the resulting minimal-rotation frame 
is invariant under rotations of the inertial coordinates 
in which the precessing waveforms are extracted, except 
for the remaining freedom of a constant overall rotation. 
Therefore, the approach proposed here results in an es- 
sentially unique adapted frame and in gravitational wave- 
forms that are similarly unique. In contrast, the imple- 
mentations of Refs. [18] and [19] choose the final rota- 
tion angle in a coordinate-dependent manner: working 
in terms of Euler angles and always setting the third Eu- 
ler angle to 7 = — by construction, this is the rotation 
about the z axis of the adapted frame. 

The key to fixing this remaining freedom by geomet- 
ric means lies in an analogy with the non-precessing bi- 
nary. In the non-precessing case, there are again many 
coordinate frames that preserve the condition that the 
radiation is primarily quadrupolar with m = ±2. One 
could arbitrarily rotate about the z axis and preserve 
this condition. In practice this does not pose a problem. 



^ More specifically, the method of Schmidt et al. calls for a rota- 
tion that maximizes the power in the (-^, m) = (2, ±2) modes. 
However, we show in Sec. II that this can be regarded — in some 
sense — as a restriction of the method of O'Shaughnessy et al. 



because it is taken as physically obvious that one would 
not analyze the waveform in a coordinate frame that is 
rotating about the z axis relative to an inertial frame. A 
fixed overall rotation about the z axis is allowed, lead- 
ing to the well-known overall freedom in the waveform 
phase. However a time-dependent rotation about the z 
axis, which could cause arbitrary frequency modulation 
in the waveform, is rejected as unnatural. 

Any frame that tracks the radiation axis of a precess- 
ing binary, on the other hand, is necessarily changing. 
As far as possible, we would like to carry over the non- 
rotating condition from the non-precessing system. In 
this case, we can describe the rotation of the frame by 
the instantaneous rotation vector w. Relative to an in- 
ertial frame, the time derivative of any vector stationary 
in the rotating frame is given by 

V = W X V . (1) 

If we denote the radiation axis by a unit vector a, we see 
that a = w X a. Taking the cross product of both sides 
of this equation with a and using the standard triple- 
product formula, we have 

a X a = {a • a) w — {a • w) a . (2) 

Using the fact that a is unit, we can rearrange this as 

w = axa-\-{a-w)a. (3) 

Of course, the component of w along a is completely un- 
determined by this equation; we need some other condi- 
tion to fix it. Now, when a = 0, as in the non-precessing 
case, we recover the natural non-rotating frame when 
a - w = 0. This is the same condition imposed by 
Buonanno, Chen, and Vallisneri [20] in the context of 
post-Newtonian template waveforms (discussed further 
in Sec. VA below). We stress the importance of this 
condition more broadly — and particularly in the context 
of numerical relativity. 

Here, we augment the methods of Schmidt et al. and 
O'Shaughnessy et al. with the condition that the instan- 
taneous rotation of the frame satisfy 

w-a = . (4) 

Hereafter, we refer to this as the condition of minimal 
rotation, as this implies that has the smallest possible 
magnitude, out of the infinitely many rotation vectors 
consistent with the known motion of the radiation axis. 
It is significant that this condition on is geometrically 
meaningful, because a is — at any instant — independent 
of the orientation of the frame in which it is found. As 
we will demonstrate below, the waveform decomposed in 
such a frame is independent of an overall rotation, up to 
a constant phase. 

Given a rotation R(t) that takes the z axis into the 
radiation axis {a{t) = R(t) z), we can use Eq. (4) to find 
a condition on R(t) that holds only if it is a minimal 
rotation. Alternatively, given any rotation that takes the 
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z axis into the radiation axis, we can easily construct 
another rotation that does the same while also satisfying 
the minimal-rotation condition. These relations are the 
key results of this paper. 

The remainder of this paper is structured as follows: 
In Sec. II we summarize the algorithms for finding the 
radiation axis presented by Schmidt et al. [18] and by 
O'Shaughnessy et al. [19]. We show that the first method 
is essentially a restriction of the second, but point out 
that with slight improvements to the numerical tech- 
niques both can be used find the correct radiation axis 
to very high accuracy — at least for simple toy models in 
which the correct axis is known. In Sec. Ill we translate 
the minimal-rotation condition, Eq. (4), into a condition 
on the rotation operator itself, and construct a method 
for imposing this condition while leaving the radiation 
axis fixed. In Sec. IV, we compare the original algo- 
rithms to this coordinate-independent method. First, we 
show that the motion of the coordinate axes is essen- 
tially invariant for our implementation, while the axis 
motion for the original methods depends sensitively on 
the orientation of the inert ial coordinate frame. We then 
demonstrate that this dependence shows up in phase the 
waveform modes using a simple post-Newtonian model. 
In Sec. V, we exhibit applications of the minimal-rotation 
frame, which demonstrate that the usual machinery used 
for non-precessing waveforms can be directly carried over 
to precessing systems in this frame. In particular, we dis- 
cuss a framework for calculating post-Newtonian wave- 
forms taking advantage of this simple frame, first pro- 
posed by Buonanno, Chen, and Vallisneri [20]. We then 
describe how to compare waveforms and construct hy- 
brids. Finally, in Sec. VI, we close with discussion of the 
benefits of this method and potential applications in an- 
alytic constructions. Two appendices detail our conven- 
tions, list some crucial formulas for rotations, and repeat 
our main results in the language of quaternions. 



II. LOCATING THE RADIATION AXIS 



to do this have been presented in the literature. We now 
review these, showing that they can be expressed in very 
similar ways, noting that both can be implemented nu- 
merically to achieve very high accuracy, and highlighting 
the crucial degeneracy present in both. 



A. The two methods 

The first algorithm, presented in Schmidt et al. [18], 
finds a frame in which the amplitudes of the (£, m) = 
(2, ±2) modes of the gravitational-wave field are maxi- 
mized. Reference [18] uses the Newman-Penrose Weyl 
scalar ^4, though the principle is the same for any radi- 
ation field — for example, the metric perturbation h. We 
will use the generic symbol q to establish the method. We 
regard q as the quantity measured in an inert ial frame, 
and q the function under a specified rotation.'^ Each of 
these can be decomposed in spin-weighted spherical har- 
monics [Eq. (A 16)], with weights of the modes denoted 
qi,m (J^'^. The relation between q^''^ and g^'^ is 
given by Eq. (A17). 

The basic idea is to find a rotation R(a,/3,7) to max- 
imize the quantity^ 



m=±2 



= E 



m=±2 



(2) 

where the V^i ^ are given in terms of the Euler angles 
by Eq. (A12). By considering the relationship between 
the coordinate systems, Eq. (A15), we can see that the 
radiation axis is given by a = R(a, /3, 7) z. 

O'Shaughnessy et al. [19] introduced another method, 
which finds an axis associated with the quadrupolar part 
of the radiation field. They begin by defining 



(5) 



The gravitational waves radiated from a compact bi- 
nary are typically decomposed in a spin- weighted spher- 
ical harmonic expansion of the field on a sphere. For 
a binary with orbital angular momentum L along the 
z axis, the dominant modes in this expansion are the 

m) = (2, ±2) modes. When L is not along the z axis, 
however, the various modes will mix, and other modes of 
the i = 2 component can dominate. For precessing bina- 
ries, this misalignment of the angular momentum and the 
z axis of an inert ial frame is inevitable, complicating com- 
parisons between simulations produced with even slightly 
different initial conditions. Moreover, the amplitude and 
phase of the modes themselves will become rapidly vary- 
ing functions of time, complicating analysis of the wave- 
forms. Both of these complications can be eliminated by 
decomposing the waveform in a non-inertial frame that 
somehow tracks the motion of the binary. Two methods 



(6) 



where La is the usual angular-momentum operator [21], 
and for simplicity of presentation we set J \qf = 1. 



^ To be precise, we define q to be the function satisfying 
q(i9',(p') = q(i9,(^), for any angles related by Il{0,i9' , cp') = 
R(q;, /3, 7) R(0, 1?, ip), where each R is a rotation operator param- 
eterized by the Euler angles as described in Appendix A. These 
conventions affect details of later results — for example, the form 
of Eq. (5) and its independence of 7 (as opposed to a). 

4 por general systems in general orientations, relations like the 
usual q^'-^ = need not hold. As a result, both terms 

in the sum over m are required, to avoid mixing of the \m\ ^ 2 
modes. 
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The radiation axis is then defined as the dominant prin- 
cipal axis of this matrix — the eigenvector with the eigen- 
value of largest magnitude. This problem can be solved 
directly with standard algebraic techniques. 

We find it useful to think of this method in a sec- 
ond way. Basic results from linear algebra show us that 
there exists a rotation operator R such that the matrix 
R R~^ is diagonal, and that the final column of 

this diagonalized matrix is the dominant principal axis. 
To put it another way, then, this method can be regarded 
as finding a rotation operator that maximizes the z-z 
component of the rotated matrix, in which case the radi- 
ation axis is just a = Rz. 

The similarity between the two methods becomes clear 
when we expand this rotated matrix and take the z-z 
component: 

(R (L(^L5))R-^)^^ 

£,m,m' 

£,771 

O'Shaughnessy et al. suggest the possibility of limiting 
the range of the sum to just i = 2. If we further limit the 
sum to m = ±2, we have = 4 times the quantity Q 
given in Eq. (5), and the method of [19] reduces to that 
of [18]. 

Important differences remain between the implemen- 
tations possible with the two methods, however. When 
O'Shaughnessy et al. sum over all relevant m modes, they 
are rotating the £ components of the waveform, which are 
geometrically meaningful. Thus, the full matrix {L(^aLb)) 
as they define it is a tensor and therefore obeys stan- 
dard rotation rules. If we limit the sum over m modes, 
we have a quantity that does not behave properly under 
rotations. This difference means that R can be solved 
for algebraically in the method of O'Shaughnessy et al.^ 
while the method of Schmidt et al. requires a more active 
maximization procedure. 

Nonetheless, we note that the method of Schmidt et al.^ 
if implemented carefully, can be made quite accurate and 
efficient. Because the right-hand side of Eq. (5) (and even 
its derivatives) can be easily expressed as a known ana- 
lytic function of the angles a, /3, 7, the problem is per- 
fectly suited to numerical optimization. We find that it is 
very easy to implement, with the code converging to the 
correct radiation axis within roughly 10~^ rad, typically 
using fewer than 10 function evaluations. This method is 
also quite robust, requiring no initial guess for the radia- 
tion axis. The speed and accuracy of this code, then, are 
essentially the same as the speed and accuracy of code 
implementing the method of O'Shaughnessy et al. 



B. Degeneracies 

In the discussion above, we glossed over a pair of de- 
generacies present in both of these methods. The first is 
trivial: the radiation axis produced by either method is 
really a directionless axis, rather than the directed axis 
we have assumed. Roughly speaking, this means that a 
may be either parallel or anti-parallel to the orbital angu- 
lar momentum. This degeneracy may be resolved by any 
convenient means, such as comparison with the coordi- 
nate angular velocity. In the following, we assume that a 
is chosen to lie parallel to the angular momentum or — at 
least — points in a consistent direction from moment to 
moment. 

The second degeneracy, however, exhibits a significant 
flaw in the methods as presented: both are invariant un- 
der rotations about the radiation axis, and therefore do 
not fix the rotation w{t) uniquely. We can see this ex- 
plicitly by looking at the behavior of the modes under 
such a rotation: 

^^i,m^imj , (8) 

where 7 is the angle of the rotation. Using this in either 
Eq. (5) or Eq. (6), we see that the phase factor cancels 
out, leaving no change to the expressions. We need to 
impose another condition to make these problems well 
posed. Both Schmidt et al. and O'Shaughnessy et al. 
break the degeneracy by simply setting the final Euler 
angle of the rotation to 0. In our conventions, this means 
setting 7 = at all times. But this choice means that 
the rotation R(a, /3, 7) depends on the inertial frame with 
respect to which the Euler angles are defined. For general 
precessing systems, it will affect the phase of the final 
waveform in highly nontrivial ways, as we demonstrate 
in Sec. IV B. 

Nonetheless, the choice of 7 = does make the par- 
ticular problem of finding the radiation axis well posed. 
In the next section we take that radiation axis, and use 
the freedom in 7 to construct a geometrically meaningful 
frame. This requires abandoning locality in time: while 
the methods of Refs. [18] and [19] can be applied for each 
t separately, our method will result in an ordinary differ- 
ential equation for the rotation matrix. 

III. MINIMIZING ROTATION 

The techniques just described give us one particular 
rotation Rra(^) = that aligns the inertial 

frame with the radiation axis. But the previously noted 
freedom in 7 means that we can first perform a rotation 
R^(t) by an angle 7 about the z axis without affecting 
the radiation axis. We now construct a rotation 

R(t) =Rra(t)R7(t) , (9) 

and solve for R^(t) such that the minimal-rotation con- 
dition, Eq. (4), is satisfied. The new rotation R(t) will si- 
multaneously satisfy the minimal-rotation condition and 
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align the inert ial z axis with the radiation axis. To find 
R^, we express Eq. (4) in terms of the rotation opera- 
tor R alone, making use of generators in the Lie algebra 
50 (3), and various relations noted in Appendix A. We 
then apply this to the case where R is decomposed as in 
Eq. (9), allowing us to solve for the minimal rotation. 



A. The minimal- rot at ion condition in terms of the 
rotation operator 

We begin by defining the equivalents of the instanta- 
neous rotation axis w and the radiation axis a using the 
isomorphism a which maps 3- vectors into 50 (3), given 
in a Cartesian basis by Eq. (A2). We write 77 = <j(t^) 
and A = a {a). Now, the dot product can also be de- 
fined for elements of 50 (3) [as —1/2 times the trace of 
the product matrix; see Eq. (A4)], allowing us to rewrite 
the minimal-rotation condition as 

n-A = 0. (10) 

Here, IJ is unknown, and A is time dependent. Therefore, 
we now translate these into expressions in terms of the 
rotation operator and the basis element Z = (j{z). 

The formula for A is simple. Recall from Sec. II that 
a = Hz. In terms of generators, this is A = RZR~^. 
The formula for 77 can be found by considering Eq. (1), 
applied to any vector v that is stationary in the rotat- 
ing frame. If we define vq := R~^(0) 'iT(O), we can write 
v{t) = Il{t)vo in the inertial frame. In 50 (3), this is 
written V = RyoR~^- Then Eq. (1) becomes 

^ (Rl/oR-') = [n.KVoK-'] (11) 
= [RR-\RVbR"^] , (12) 

where the second line comes from expanding the deriva- 
tive using Eq. (A6). If equality is to hold for arbitrary 
we must have 



B. Solving for the initial rotation 

To find R satisfying the minimal-rotation condition, 
we now insert Eq. (9) into Eq. (15). Assuming that Rra 
is known, and using R^ = exp(7Z), this will give us a 
condition on 7. First, we calculate 

R-^R = e-^^R-^Rrae^^+7Z . (16) 

Note that conjugation by ^ does not affect the com- 
ponent along Z. Therefore, plugging this result into 
Eq. (15) and rearranging, we obtain 

7=(-R-'Rra)-^. (17) 

Because Rra is known, we can simply evaluate the right- 
hand side, integrate to find 7(t), and insert this back into 
Eq. (9) to find a minimal-rotation operator that takes the 
z axis into the radiation axis. 

We emphasize that the derivations of Eqs. (15) 
and (17) did not assume any features of R and Rra other 
than the fact that they rotate the z axis of the iner- 
tial frame into the radiation axis. In particular, we did 
not assume that the final Euler angle was zero — or in- 
deed use any expression in terms of the Euler angles. 
Also, though inspired by the Euler angles, the defini- 
tion R^ := exp(7 Z) is independent of coordinates on 
SO (3).^ Therefore, Eq. (17) is a geometrically general 
equation: it should hold regardless of any coordinates 
we might choose for 5*0(3), and should take the same 
form regardless of the inertial frame we use to find the 
radiation axis. One result of this is the fact that frames 
satisfying the minimal-rotation condition are unique up 
to a constant overall rotation about the radiation axis, 
corresponding to the integration constant obtained from 
integrating Eq. (17). We discuss this freedom further in 
Sec. VB. 

Nonetheless, for the purposes of implementation, an 
explicit formula involving the Euler angles will be use- 
ful. When Rra(^) = R(a(t), /3(t), O) , a straightforward 
calculation using Eqs. (A7) and (A8) gives us^ 

7 = —a cos f3 . (18) 



i7 = RR-^. (13) 

Now, using these expressions for A and 77 in Eq. (10) 
and rearranging a little, we get another form of the 
minimal-rotation condition: 

(RR-^) • (RZR-i) =0 . (14) 

This is precisely the minimal-rotation condition of Eq. (4) 
in operator form. We can simplify this expression 
slightly. Noting that the dot product is invariant un- 
der rotations, we apply the inverse rotation to each part 
of the product, obtaining 

(R-^R) 'Z = . (15) 



^ Note the distinction between the choice of basis for and the 
choice of coordinates for SO(3). Here, the basis (x,y,z) pro- 
duces the canonical basis of generators (X,Y,Z) for 50 (3) by 
eigenvector problems. In particular, Z represents the unique 
generator having eigenvector z with eigenvalue 0, and eigenvec- 
tor {x-\-iy)/\/2 with eigenvalue i. Also, the exponential function 
is defined geometrically (see Eq. (Al) or Ref. [22]), which shows 
that exp(7Z) is independent of coordinates on SO(3). Simi- 
larly, the arbitrary choice of z axis will not affect the form of 
these equations; we can choose any unit vector as the z axis, 
and as long as the rotation operators take that vector into the 
radiation axis, the expressions will not change. 

^ Near the coordinate singularities at /5 = and /3 = tt, the 
value of a will contain substantial numerical noise. Differ- 
entiating a, as in this equation, simply magnifies that noise. 
For some configurations, we find improved numerical results 
when integrating this equation by parts and implementing it as 
7 = —a cos /3 — / cx $ sin /3 dt. 
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Schmidt et al [18] pointed out that the orbital frequency 
in the rotating frame given by Rra should be roughly 
Q -\- a cos P (in our notation), where Q is the magnitude 
of the orbital frequency measured in the non-rotating 
frame. Thus our adjustment to their technique can be 
thought of as removing that second term, so that the or- 
bital frequency in the rotating frame given by R should 
be roughly Q. We will return to this observation in 
Sec. V A to discuss a simplification of the post-Newtonian 
representation of precessing systems. 



IV. THE EFFECTS OF ENFORCING THE 
MINIMAL-ROTATION CONDITION 

No preferred inertial frame exists for general precess- 
ing systems. Indeed, no preferred axis exists for choosing 
an inertial frame. When the total angular momentum J 
points in a constant direction, this can be a useful choice 
of axis, but J will change direction for inspir ailing pre- 
cessing systems. We might therefore expect data from 
different numerical simulations, for example, to be pre- 
sented in different inertial frames. A key concern, then, 
is the behavior of the rotated frame under fixed rota- 
tions of the inertial frame. We now demonstrate that a 
frame aligned with the radiation axis using 7 = behaves 
poorly under such rotations, whereas such a frame with 
the minimal-rotation condition imposed is essentially in- 
variant. First, we observe the motion of the rotated 
axes in the two cases. Then, we inspect the behavior 
of the phase of the waveform decomposed in the rotating 
frames. 



A. Rotation history in different inertial frames 

For any rotation R(t), we visualize its "rotation his- 
tory" by plotting the paths of the tips of the rotated basis 
vectors, x', and i', on the unit sphere. We demon- 
strate the rotation histories for a toy model in which the 
radiation axis a precesses about a fixed axis F, where a 
is inclined to F by an angle of 25°. In the first panel 
of Fig. 1, we show the rotation history for this system 
when F is along the z axis and the final Euler angle is 
simply set to 7 = 0. The blue circle traces the path 
of the radiation axis — which coincides with the basis 
vector — while the red and green curves trace the paths 
of the x' and y' axes, respectively. 

If, instead of being aligned with the z axis of the iner- 
tial frame, the F axis is tipped, we obtain different rota- 
tion histories. Later panels of Fig. 1 show the histories 
when F is tilted by the given amount. In each panel, the 
blue curve remains the same, being simply shifted on the 
sphere by the given tilt. However, the paths of the other 
two axes of the adapted frame change drastically as the 
inclination of F is changed, even undergoing topological 
transitions as this inclination passes the 25° precession 



tilt: 



tilt: 10° tilt: 20° tilt: 30° 
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tilt: 40 ° tilt: 50 ° tilt: 60 ° tilt: 70 ° 

tilt: 80° tilt: 90° 



FIG. 1. Rotation histories for a simple precessing sys- 
tem when 7 = 0. The axis denoted by the blue curve (the 
small, perfectly circular loop) precesses through a single cycle 
around a fixed axis that is tilted by varying amounts. The 
other axes of the tracked frame (red and green curves) take 
paths that depend on the inclination of the axis around which 
the precession occurs, relative to the inertial frame. The y 
axis, in particular, is forced to remain on the x-y plane of the 
visualization coordinates, by the choice of 7 = 0. If we had 
defined Euler angles by z-x-z rotations, as in Ref. [18], then 
it would have been the x axis that was forced to remain on 
this plane. 



inclination and its 65° complement. The tracked frame 
thus has a time-dependent rotation about the radiation 
axis, which will show up as a time-dependent modulation 
of the waveform phase. This modulation is determined 
by the choice of inertial frame, and therefore the phase of 
a waveform measured in a frame obtained with 7 = is 
not invariant in any useful sense. This phase modulation 
will be examined directly in Sec. IV B. 

We can repeat this comparison of rotation histories 
when the third Euler angle is set by the minimal-rotation 
condition, Eq. (18). The results are shown in Fig. 2. In 
this case, the rotation histories have the same shapes, but 
are simply tilted with respect to each other. That sim- 
ilarity shows that frames constructed with the minimal- 
rotation condition are essentially^ invariant under fixed 
rotations of the inertial frame. As we will now see, the 
waveform measured in such a frame is similarly invariant. 



^ There is still an overall freedom in a fixed rotation about the ini- 
tial radiation axis. This corresponds to the standard ambiguity 
in orientation, which must be fixed by other methods, discussed 
in Sec. VB. 
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FIG. 2. Rotation histories for a simple precessing system 
when 7 is set by Eq. (18). These systems are just the same 
as in Fig. 1, except that the frame is given a final rotation 
about the z axis to satisfy the minimal-rotation condition. 
The paths are identical in each case, but tilted with respect 
to each other. This shows that the minimally rotated frame 
is essentially invariant under fixed rotations of the inertial 
frame. Note that the paths in these figures represent many 
precession cycles, each corresponding to a single cycle of the 
"zigzag" pattern. 



B. Waveforms in different inertial frames 

The methods set forth in Refs. [18] and [19] fix the 
axis of the rotated frame in an invariant way, but rely on 
coordinate-degrees of freedom to fix rotations about that 
axis. Such a rotation translates directly into the phase of 
the waveform measured in the rotated frame. Since the 
minimal-rotation condition is imposed with a rotation 
about this radiation axis, only the phase of the waveform 
is affected. Therefore, we ignore the waveform amplitude, 
as it is invariant as soon as the radiation axis is fixed. We 
examine the waveform phase with two examples — a sim- 
ple analytical example first, followed by a more realistic 
model. 

Imagine a system with very small precession, where the 
radiation axis moves with constant speed along a vanish- 
ingly narrow cone centered about the z axis. Because the 
cone is so narrow, we might hope that the phase of the 
waveform q measured in the inertial frame will be close 
to the phase measured in the rotated frame. Assuming 
the radiation axis is in the x-z plane at time t = 0, the 
rotation found by the methods of Refs. [18] and [19] will 
be R(at,/3,0), for constants a and p. We can relate the 
modes of the waveform measured in the rotated frame 



to the modes measured in the inertial frame q ' 



using Eq. (A17b). If we approximate /3 ~ 0, then the 

Wigner matrices are nonzero only for m = m\ 

and we have 



That is, the waveform acquires an additional linearly 
increasing phase in the rotated frame. If the complex 
phase of this mode is the change in going from a 

non-precessing system to a system precessing with a very 
small opening angle is 



(19b) 



Note that the additional phase only depends on the num- 
ber of times the system has precessed, not the size of the 
precession angle. If, on the other hand, we impose the 
minimal-rotation condition, Eq. (18) gives us 7 —at, 
which cancels the additional phase, so that g^'^ ~ q^''^. 
In this sense, the minimal-rotation frame is much more 
natural. 

More importantly, however, a frame with minimal ro- 
tation behaves nicely under fixed rotations of the inertial 
frame. For example, we take the same system as above, 
but tilt the inertial frame slightly so that the precession 
cone lies close to, but does not contain the z axis. In 
this case, the rotation will be Il{(p, 0), where (1?, cp) are 
the usual spherical coordinates of the axis, which we can 
calculate by simple trigonometry. Using this result to 
transform the modes, we find 



^£,m ^ £,171—1171 arccot[cot(Q! csc(q! t)] 



(20a) 



where a is the same constant as above and > 1 is 
a constant that depends on the particular values of the 
precession and the tilt angles.^ In this case, the change 
to the waveform phase is bounded, but oscillatory: 



/'^(t) ^ 0''^(t) - m arccot[cot(dt) + 6>csc(dt)] . 

(20b) 

Thus, a slight change in the inertial frame causes a dras- 
tic change in the behavior of the waveform modes, which 
is associated with the topological change seen in Fig. 1 
when the tilt exceeds 25°. In the minimal-rotation frame, 
on the other hand, the phase change in Eq. (20a) is 
counteracted by the adjustment to 7, and we still have 

These features also show up in systems with significant 
precession. We now turn to a post-Newtonian model of 
such a binary. The system we choose has equal-mass 
black holes, with spins x = 0-99 initially parallel to each 
other and orthogonal to the orbital angular momentum. 



, £,m — i mat 
q e 



(19a) 



In fact, this form of the equation applies for all small precession 
angles and tilts. If B is the tilt angle, then 6 = Bf/S, so = 
corresponds to the case where the precession cone is centered on 
the z axis, and Eq. (20) reduces to Eq. (19). Similarly, < ^ < 1 
corresponds to a small tilt, for which the z axis is still within 
the precession cone. In this case, the additional phase is roughly 
linear, with a superposed oscillation. = 1 corresponds to the 
case where the radiation axis passes through the z axis, which is 
the coordinate singularity of the Euler angles, meaning that the 
waveform phase in the 7 = frame is actually undefined. Finally, 
> 1 corresponds to a tilt that is larger than the precession 
angle, so that the z axis is not enclosed in the precession cone. 
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FIG. 3. Waveform frequency in various frames aligned with 
the radiation axis. Here, waveform frequency refers to the 
time derivative of the phase of the {i,m) = (2,2) mode of 
the gravitational waveform. The solid red line shows the fre- 
quency measured in a frame derived from the inertial frame 
by a rotation in which the third Euler angle 7 is set to 0, 
while the solid blue line shows the same quantity in a frame 
for which 7 satisfies the minimal-rotation condition. Clearly, 
the latter curve is much smoother. We also show as dotted 
lines the same quantities when the physical system is tilted by 
10°. The dotted blue line coincides with the solid blue line, 
showing the invariance of the waveform in that frame. 



FIG. 4. Change of phase measured in frames aligned to the 
radiation axis when the physical system is tilted by 10°. This 
phase difference is defined by Eq. (21), and plotted for two 
cases: the first (in red) where 0^'^ refers to phases measured 
in a frame obtained from the inertial frame with a rotation 
in which the final Euler angle 7 is set to 0; the second (in 
blue) where 0^'^ refers to phases measured in frames satis- 
fying the minimal-rotation condition. The change in phase 
for the minimal-rotation frames is to within numerical error 
(roughly a part in 10^ here). Note that this waveform extends 
for roughly the length of time such a system would be in the 
sensitive band of Advanced LIGO if the total mass were about 
IOM0. 



which initially coincides with the z axis. The precession 
cone has an opening angle of roughly 15° initially, grad- 
ually widening to about 21°. Using the post-Newtonian 
waveform, we can the find the radiation axis with the 
methods described in Sec. II, then decompose the modes 
of the waveform either in a frame with 7 = or in a 
frame with minimal rotation. 

Again, we see the two features identified above. First, 
the waveform decomposed in the minimal-rotation frame 
appears to be smoother than the waveform decomposed 
in the 7 = frame. In particular, while the amplitudes 
are identical in the two frames, the phase of the (£, m) = 
(2,2) mode in the 7 = frame is constantly increasing 
relative to the phase in the minimal-rotation frame, and 
jumps each time the radiation axis passes near the z axis 
(each time a cos P is large). We plot the frequency of the 
(2, 2) mode measured in the two frames in Fig. 3, where 
the phase jumps show up as spikes. 

Second, the waveform phase in the minimal-rotation 
frame is invariant (up to a constant) under overall ro- 
tations of the inertial frame in which the waveform is 
measured — which is not the case for the 7 = frame. We 
illustrate this by tilting the post-Newtonian system by a 



10° rotation about the y axis, and redoing the decompo- 
sition in the two frames. The frequencies for this rotated 
system are also plotted (as dotted lines) in Fig. 3, where 
we see that the curves for the minimal-rotation frame lie 
on top of each other, while the 7 = curve changes but 
still exhibits spikes. 

Figure 4 plots the phase-differences of the (2,2) mode 
between the two post-Newtonian evolutions differing by 
a rotation by 10°, 



2,2 _ 



/.2,2 

^untilted 



^2,2 
^tilted- 



(21) 



The minimal-rotation frame is invariant under this rota- 
tion, and indeed the phases are identical to within the 
numerical error (as measured by convergence of the in- 
tegration of Eq. (18)). The coordinate-dependent choice 
7 = 0, however, results in phase differences of multiple 
gravitational-wave cycles. 

While we designed this example to be a rigorous 
test of the frame- alignment techniques, the later stages 
of inspiral and merger provide an even more stringent 
test, as precession frequency scales with orbital velocity. 
The minimal-rotation condition provides an easily imple- 
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merited solution to the problems presented by the 7 = 
frame. 



V. FURTHER APPLICATIONS OF THE 
MINIMAL-ROTATION FRAME 

The minimal-rotation frame aligned with the radiation 
axis allows us to describe gravitational waveforms from 
precessing binaries very nicely. The amplitude and phase 
are smoothly varying functions, and are invariant (up to 
an overall phase offset) under rotations of the inertial 
frame. The significance of this construction is broader 
than one might immediately imagine. We now discuss 
some of the most important uses of gravitational wave- 
forms, and show that the minimal-rotation frame can 
make the necessary manipulations much simpler. 



A. Post-Newtonian waveforms 



Assuming that Ln is given by the angles a and ^5, then 
as in Sec. IIIB we can define the rotating frame by the 
operator 

R(a(t),/3(t),-/a cos/3dt) . (24) 
In such a frame, Eq. (22) becomes 

^inertial ~ ^orbital • (^^) 

This means that the pN orbital elements can be inte- 
grated in the usual way. That is, the orbital phase ^ 
obtained by integrating the angular frequency is just the 
phase in the rotating frame. This can then be inserted 
into the standard expressions [23, 25-27] for the wave- 
form modes, ^ which gives the waveform in the rotat- 
ing frame. Together with the rotation operator given 
by Eq. (24), this describes the waveform completely. In 
particular, it can readily be transformed to the inertial 
frame using Eq. (A17). 



Buonanno, Chen, and Vallisneri [20] proposed the con- 
struction of template waveforms for precessing binaries, 
expressing the metric perturbation using the minimal- 
rotation frame (which they simply called the "precessing 
frame"). Their motivation was to separate the response 
of a detector into two parts: one describing the intrin- 
sic qualities of the waveform; the other describing the 
position and orientation of the detector. 

This description differs from the more usual ap- 
proach [15, 16, 23, 24] of tracking the binary's orbital 
parameters in a frame related to the inertial frame by a 
rotation taking the z axis into the Newtonian angular mo- 
mentum I/N, with the final Euler angle set to 7 = 0. This 
results in a nontrivial relationship between the orbital 
frequency in the inertial frame, l^inertiai, and the orbital 
frequency in the instantaneous orbital plane, l^orbitai- 

^inertial = ^orbital + ^ COS /3 . (22) 

In such a framework, this equation must be integrated as 
part of the post-Newtonian system, along with the equa- 
tions for the angular momentum, spins, and the usual 
fiux and orbital energy, resulting in the motion of the bi- 
nary as a function of time. This motion is then inserted 
into expressions for the gravitational waveform in the in- 
ertial frame. Because of the complicated motion, these 
expressions are necessarily even more elaborate than the 
standard expressions for motion in the orbital plane [23] . 

We can substantially simplify the analytical prescrip- 
tion by evolving the orbital elements in the inertial frame, 
but writing the waveform in terms of spin-weighted 
spherical harmonics in the rotating frame aligned with 
the angular momentum [15]. This is particularly easy 
if the frame satisfies the minimal-rotation condition, 
Eq. (4). In this case, the instantaneous rotation vector is 

= I/N X Ln . (23) 



B. Comparing waveforms and constructing hybrids 

Having obtained some model waveform, one of the first 
things we typically do is compare the result to some other 
waveform. For example, we might run a numerical sim- 
ulation at two different resolutions or using two different 
numerical codes, and compare the output waveforms to 
get an estimate for their accuracy. Or we might construct 
an analytical model, and compare it to a numerical model 
of the same system. And of course, numerical simulations 
can only describe the last portion of the inspiral of a bi- 
nary system, so we frequently combine analytical and 
numerical waveforms into "hybrid" waveforms. In each 
of these cases, the first step is to ensure that the two 
waveforms are expressed in the same coordinate system, 
which is generally referred to as alignment. 

With various simplifying assumptions, alignment boils 
down to setting the time offset and the relative orienta- 
tion of the frames in which the waveforms are measured. 
For a non-precessing binary, as discussed in Sec. I, this 
further reduces to setting the time offset and a single 
rotation around the z axis. That is, one waveform is 
adjusted as 

/i^'^(t) ^/i^'^(t + At)e-^^^'^ , (26) 



^ This prescription is sufficient at the level of knowledge of post- 
Newtonian waveforms given in the references. However, by treat- 
ing the waveform as if it were determined only by instantaneous 
positions and velocities, it neglects contributions which may be 
relevant at higher orders. For example, the standard expressions 
for the waveform amplitudes assume accelerations are orthogonal 
to the orbital angular velocity, which is not the case for precess- 
ing systems. Nonetheless, these contributions can be calculated 
and projected into the rotating frame. As these are higher-order 
spinning terms, we ignore them at this level. 
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for time offset At and phase offset A^. The criteria used 
to choose those offsets typicahy ensure that the phase 
and frequency of /i^'^ are the same in both waveforms at 
a particular instant, for example. 

Precession complicates this simple picture, because we 
can no longer rotate the second system by just one angle 
A^. In general, the inertial frames will be related by 
some rotation Rf, which rotates the second frame into 
the first. We would need to solve for this rotation, which 
might involve rotating the entire i = 2 component of the 
waveform to minimize some measure of the difference 
between the waveforms, for example. This is cumber- 
some, and would inherently depend on the inertial frame 
in which the waveforms are measured. 

Alternatively, we can measure the two waveforms in 
minimal-rotation frames aligned with their respective ra- 
diation axes. In that case. Sec. IV B showed that, again, 
the only freedom to choose coordinates for the waveforms 
lies in the time offset At and a single phase offset rotat- 
ing the system about the radiation axis by A^. This, of 
course, would not be the case if the frame failed to satisfy 
the minimal-rotation condition. 

These waveforms will also come with rotations Ri(t) 
and R2 {t) , which describe the orientation of their frames 
relative to some inertial frames. Again, these inertial 
frames need not be the same. However, if our criteria for 
setting the time and phase offsets make the waveforms 
equal at some fiducial time tfid, we can now directly solve 
for their relationship: 

Rf = Rr'(tfid)R2(tfid + At)e^*^ . (27) 

This value of Rf is then a constant, to be used at all 
times. More generally, we define a new rotation operator 

K'^it) = R2(^ + At) Rf"^ . (28) 

This operator is more directly comparable to the first 
rotation operator, in the sense that Ri(tfid) = R2(tfid). 

Now, waveform comparisons for non-precessing sys- 
tems typically measure differences between amplitudes 
and phases of the aligned waveforms. In precessing sys- 
tems decomposed in the aligned frame, it makes sense to 
show the same quantities. However, it is now also im- 
portant to examine how the frames are aligned as func- 
tions of time. The first interesting quantity to compare 
might be the angle between the two radiation axes, ai 
and a2 = R2 z. A simple formula that expresses this 
angle is 

:= arccos [(Riz) ■ {^2^)] . (29) 



Many other possible criteria exist for choosing these offsets (see, 
e.g., Ref. [28]) but the effect on the waveform should always be 
given by Eq. (26). Therefore we ignore the particular criteria in 
use. 



We also need to understand the relative rotation about 
the radiation axis, which is particularly important be- 
cause it translates directly into the phase of the wave- 
form. There is no unique way to define this phase, when- 
ever the two radiation axes are misaligned. In anal- 
ogy with the above, we might use arccos [(Ri y) • (R2^)]. 
However, this can be dominated by the tilt of the ra- 
diation axes, so we attempt to remove that part of the 
rotation. That is, we define the rotation Ri^2' that takes 
ai into a2 "directly" by rotating through an angle —Sf3 
about the axis ai x This can be conveniently cal- 
culated by the axis-angle formulation or quaternions, as 
described in Appendix. B. Then, the following fulfills our 
needs: 

:= arccos [(Ri^2' Ri y) ' (Rs^)] • (30) 

These two angles are enough to characterize the misalign- 
ment of two frames. 

We can also use the aligned frame to construct hybrid 
waveforms joining analytical and numerical waveforms in 
the simple manner used for non-precessing systems [29]. 
When hybridizing, we use information from only the first 
waveform before some time ti, and information from only 
the second waveform after some time ^2, with some tran- 
sition in between. This is typically accomplished using 
a transition function r{t) that equals 1 before ti, after 
^2, and transitions smoothly in between. Then, for any 
quantity such as the amplitude or phase of a particular 
mode, we define that quantity in the hybrid as 

ghybrid(t) = qi{t) T{t) + q2{t) [1 - r(t)] . (31) 

The same method does not apply trivially when applied 
to the rotation operators Ri and R2, but we can suggest 
a simple method. Various techniques have been devel- 
oped by researchers in computer graphics to allow inter- 
polation of rotation operators — typically with such unfor- 
tunate names as slerp [30], nlerp, and even log-quaternion 
lerp [31]. Regardless of the details, any of these can be 
used to define the interpolant Rinterp(^; ^), which equals 
Ri(t) when x = 1 and R2(^) when x = 0. We can then 
define the hybrid rotation operator as 

Rhybrid(^) = Rinterp('?"(t); t) , (32) 

completing the formulation of the hybrid. 

VI. DISCUSSION 

A new element has now been added to the description 
of a gravitational waveform. The complete description 
consists of three elements: 

1. specification of the inertial frame in which the 
waveform may be measured, 

2. the operator R(t) that rotates the inertial frame 
into the frame in which the modes of the waveform 
are decomposed, and 
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3. the modes of the waveform as measured in this ro- 
tated frame. 

The new element is item 2; previously, the waveform 
would simply be decomposed in the inertial frame. Ref- 
erences [18] and [19] introduced criteria for choosing R(t) 
such that the amplitudes of the waveform modes become 
simpler. Our contribution has been the introduction of 
an additional criterion which simultaneously makes the 
choice of R(t) essentially unique and makes the phases 
of the waveform modes simpler. Crucially, we developed 
the simple formula in Eq. (18) for imposing our criterion. 

This change to the description of waveforms has many 
benefits for precessing systems. First, the amplitude 
and phase of the modes will be more nearly approxi- 
mated by smooth functions. This means that less storage 
space will be needed to describe the data from numerical 
simulations, which is an increasingly important problem 
in gravitational- wave modeling [27]. Smooth functions 
are also crucial to several of the assumptions underly- 
ing extrapolation of numerical waveforms to infinite ex- 
traction radius, for example [32]. In general, we expect 
that strong spin-spin couplings will imprint the waveform 
modes with some non-smoothness, but this method does 
at least remove the inessential coordinate-dependent fea- 
tures, as illustrated vividly in Figs. 3 and 4. 

Second, analysis of the waveforms becomes simpler. 
The machinery of waveform manipulation for non- 
precessing systems is well developed. Analytical con- 
structions, methods for comparing waveforms to demon- 
strate convergence or to measure differences from analyt- 
ical systems, and hybridization techniques are all broadly 
understood and applied. The basic approaches, how- 
ever, are not designed for precessing systems analyzed 
in an inertial frame; many of the simplifying assump- 
tions break down. The strength of this reformulation 
to include the rotation operator is that it simplifies the 
waveform modes, so that we can again use techniques 
designed for the non-precessing case. 
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Appendix A: Conventions and essential formulas 

We begin by briefiy reviewing the formalism of rota- 
tions of about the origin. Such rotations form a group. 



described by orthogonal 3x3 matrices with determinant 
+1. The group is denoted 5*0 (3), and the group oper- 
ation is composition of rotations. SO (3) also forms a 
manifold satisfying certain consistency properties — it is 
a Lie group [22]. The tangent space to that manifold at 
the point corresponding to the group's identity element 
is called the Lie algebra 50 (3), consisting of all antisym- 
metric 3x3 matrices. This Lie algebra is familiar as the 
generators of rotations: any R G 50 (3) gives rise to a 
rotation R G 5*0 (3) via 



exp : 



k=0 



(Al) 



We also have an isomorphism between 50 (3) (with the 
matrix commutator as product) and the standard 3- 
vectors (with the cross product). This map is most 
easily represented in a standard Cartesian frame, where 
the vector v = X(^]^^ is mapped as 



-v^ yy 

—v^ 

-yy 



(A2) 



In particular, we will use the isomorphism relation 

a{v X w) = [a{v),a{w)] , (A3) 

for any iT, w G V^. We win also use the Cartesian basis 
{x^y^z) for V^, corresponding to the basis (X, F, Z) for 
50 (3). We can decompose any element of 50 (3) in this 
basis by translating the dot product into 50 (3): 



(A4) 



which satisfies a -b = a{a) • cr{b). This can be useful as 
in Eq. (15), for example, where the component along Z 
must be taken. 

In this context, several formulas will be very useful. 
For any rotation R G SO (3) and vector iT G V^, 



For any differentiable curve R(t) G 50(3), 
d. 



dt 



(A5) 



(A6) 



If R{t) G 50 (3) is a curve such that R(t) = exp R{t)^ then 
whenever R{t) 7^ we can calculate 



R R — R 



1 — cos r 



[R,R] + !-^^[R,[R,R]] , 

(A7) 

where r is the magnitude of the nonzero eigenvalues of R, 
which also equals the vector norm |?^ = |cr~^(i?)|. For 
any t such that R{t) = 0, only the first term remains. 
Finally, for A, B e 50 (3) with A 7^ 0, 



A ^ A ^ 1 — COS a r . r . sin a , , 
e^Be-^ = B^ [A, [A, B]] + [A, B] 



(A8) 
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where a is similarly the magnitude of the nonzero eigen- 
values of A. Obviously, when A = 0, only the first term 
remains. 

Euler angles form one convenient set of coordinates 
for 5*0 (3). We define the Euler angles using the z-y^- 
z" convention, where the first rotation is through an 
angle a about the z axis, the second through /3 about the 
(new) y' axis, and the third through 7 about the (new) 
z" axis. Note that this is equivalent to rotations in the 
opposite order about the fixed set of axes z-y-z — which is 
an especially useful form for calculations. In particular, 
we can express the rotation operator as 



R(a,/3,7) =e"^e^^e^^ 



(A9) 



Note that, in this convention, the singularities of the Eu- 
ler angle coordinates occur when /3 = or /3 = tt. 

The Wigner matrices V^^^ form a representation of the 
rotation group. That is, we know that a composition of 
rotations in S'0(3) given by 

R(a, 7) = R(a', fi'.i) R(a", /?", 7") (AlO) 
implies the relation 

E^™'.™"(«''/5',y)P5,™(a",/?",7") . (All) 

m" 

We can find an explicit formula for the V^^^ as in Ref. [34] 
which, in our conventions, gives 



/ (£ + m)!(£- 


m)! 


(^ + mO!(^- 





gi(maW7) ^(-l)'^ ^^ + m'^ ^ 



p — m — m! 



X sm 



^ \ 2£— 2p+m+m'' f 2p—m—m' 



(A12) 



Goldberg et al. [21] showed that the spin- weighted 
spherical harmonics [35] are simply special cases of this 
expression. Adopting conventions to agree with Ref. [27], 
we have 



^J-tlvl^]_M^,^) . (A13) 

4 TT ' 



Combining Eqs. (All) and (A 13), we can immediately 
obtain the transformation law for the spherical harmonics 
under rotation of the coordinates [36]: 

,r,,^.(^^(^0 = E^m^m(«./^.7).>^^,m(^,^) , (AM) 



Note that the z-x'-z" convention is more standard in classical 
mechanics — as in Ref. [33], for example. The z-y'-z" convention, 
however, is more standard when using spherical harmonics, as 
we do here. 



for angles satisfying 

R(0,l?^(^0 = RK/^,7)R(0,l?,(^) . 



(A15) 



Note that the rotation R(a,/3,7) is uniquely fixed by 
this condition in a way that it would not be by simply 
requiring the point defined by {'d'^cp') to rotate into the 
point defined by {i^^(p). The fact that such a condition 
on the coordinates would not uniquely define R(a, /3, 7) 
is the central problem addressed by this paper. 
We decompose a function q of spin weight s as 



(A16) 



We then define the rotated function q satisfying 
q{'d'^(p') = q{'d^(p)^ where the angles are again related 
by Eq. (A15). We can decompose q into spin- weighted 
spherical harmonics in (i^^, (/?'), and use Eq. (A 14) to find 
the relations between the components: 

9'''" = E9'''"'^-U("'/5'7), (A17a) 

or equivalent ly 

qi.^ = E q''"^' -13, -a) . (A17b) 

These relations are special cases of more general trans- 
formations derived by Gualtieri et al. [17]. 

Appendix B: Main results in quaternion form 

Another useful parameterization of the rotation group 
is given by the set of unit quaternions, which are useful 
as an efficient numerical method for representing rota- 
tions. Here, we repeat the main results of this paper 
(the minimal-rotation condition and a method for im- 
posing it) in the form of quaternions. We also express 
the Wigner matrices P^^^ as functions of a quaternion, 
instead of Euler angles. 

Quaternions may be thought of in many ways, but the 
form we find convenient here is that of a scalar plus a 
vector. The notation we use will be Q = <7o + ^ = 
(go, ^2, <73), where qo is the scalar part, q is the vector 
part, and (^1,^2, ^3) are the Cartesian components of the 
vector. The conjugate of the quaternion is Q = qo — q. 
The distinctive feature of quaternions is their unusual 
multiplication rule: 

PQ = {poqo-p-q)^{poq^qop^px q) • (Bi) 

Note that this product is associative (unlike the vec- 
tor cross product) but is neither commutative nor anti- 
commutative for general quaternions. The unit quater- 
nions satisfy Q Q = 1, and can be thought of as points on 
the unit 3-sphere. The set of all unit quaternions forms 
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a group locally isomorphic to SO (3), where rotation of a 
vector V by the unit quaternion Q can be expressed by 

v' = QvQ. (B2) 

In fact, the quaternions provide a double cover of 5*0 (3), 
as can be seen above by the fact that Q and —Q produce 
the same rotation, but this does not cause any practical 
problems. The relationship with the axis- angle formu- 
lation of rotations is particularly clean: for a rotation 
through an angle 5 about an axis the quaternion is 
given by Q = cos((^/2) +sin(J/2) w. In any case, we may 
think of a general rotation R as being precisely equiv- 
alent to some unit quaternion. For further details, we 
refer to standard texts (e.g., Ref. [37]). 

A calculation very similar to the derivation of Eq. (13) 
shows us that the instantaneous rotation vector associ- 
ated to the unit quaternion R is given by 

w = 2RR. (B3) 

If R rotates the z axis into the radiation axis, a = Rz R, 
then the minimal-rotation condition given by Eq. (4) can 
be rewritten as 

{RzR)o=0, (B4) 

where the subscript on the left-hand side takes the scalar 
part of the expression. This equation is equivalent to the 
condition on the rotation operator itself given in Eq. (15). 



As before, we may impose this condition by applying 
an initial rotation about the z axis. Thus, if R^a is any 
quaternion rotation that takes the z axis into the radia- 
tion axis, we define the rotation R = R^^R-f^ where Rj 
is a rotation through an angle 7 about the z axis. Then 
R satisfies the minimal-rotation condition if 

7 = 2(iira^^ra)0 • (B5) 

Again, given the result R^a of some axis-alignment 
method, we can evaluate the right-hand side, integrate, 
and construct the total rotation. 

Finally, if the result is to be used to rotate modes of 
a waveform, we need to express the Wigner matrices in 
quaternion form. This is done most simply by relating 
the Euler angles to various components of the quaternion, 
and substituting the results in Eq. (A12). We find 

T^W (r>^(_n^m-i I (^ + m)!(£ 

^m',mW I ^) ^ (i ^ rny.{i - my. 

X [{Rz)s - i {Rz)or^^' [{Rz), + i {Rzhr-^^ 

p \ P J \p — m — m' J 

X [l-{RzR),]'^'^'-'' [l + iRzR),]"—' . (B6) 

This expression can be used in Eq. (A17a); for 
Eq. (A17b), R = R~^ should replace R in this expression. 
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